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I Abstract. We show that the soft X-ray spectra and light 

■ curves observed with the ROSAT and EUVE from the 
\ closest known millisecond pulsar J0437-4715 can be in- 

■ terpreted as thermal radiation from two hot polar caps 
' whose emitting layers (atmospheres) are comprised of hy- 

■ drogen. The simplest model yields a uniform temperature 
; of (0.8 - 0.9) X 10^ K within a cap radius of 0.7 - 0.9 km. 

■ The spectral fits indicate that the temperature may be 
nonuniformly distributed along the cap surface. The dis- 
tribution can be approximated by a central core heated 

\ up to (1 — 2) X 10® K within a radius of 0.2 — 0.4 km, sur- 
. rounded by a colder rim with temperatures (3 — 5) x 10^ K 
' extending out to 2— 6 km. The polar cap interpretation im- 
plies low column densities, (1 — 3) x 10^^ cm~^, and a high 
degree of ionization, > 20%, of the interstellar hydrogen 
towards the pulsar. The inferred bolometric luminosity of 
' the polar caps, (1.0 — 1.6) x 10^° erg s~^, is in excellent 
agreement with the predictions of the slot-gap model of ra- 
dio pulsars developed by Arons and his coworkers. Similar 
. polar cap radiation should be emitted by other millisecond 
'■ pulsars, although in some of them (e. g., PSR B1821-24) 
the soft X-ray flux is dominated by the nonthermal radi- 
ation from pulsar magnetospheres. 

' Key words: stars: neutron - X-ray: stars - pulsars: indi- 
vidual: PSR J0437-4715 



1. Introduction 

The bright 5.75 ms pulsar J0437-4715 was discovered by 
Johnston et al. (1993) during the Parkes survey for mil- 
lisecond pulsars. The pulsar is in a 5.74 d binary or- 
bit with a cool (Tcoior 4000 K), low-mass (~ 0.2Mq) 
white dwarf companion and is surrounded by a bow-shock 
nebula (Becker et al. 1993; Bailyn 1993; BeU, Bailes & 
Bessel 1993; Danziger, Baade & Delia Valle 1993; Bell 
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et al. 1995). It is an old object, with a characteristic age 
T = P/2P ~ 5 X 10^ yr, low magnetic field B - 3 x 10* G, 
and rotational energy loss E = 4 x 10^^^ erg s^^. The ob- 
served dispersion measure of 2.65 pc cm~^ implies a dis- 
tance d ~ 100 — 180 pc, making this the closest known 
millisecond pulsar. Sandhu et al. (1997) reported the dis- 
tance d ~ 178 ± 26 pc, from parallax measurements. The 
pulsar shows significant radio emission over at least 80% 
of the pulse period, with a complicated mean pulse shape 
varying with radio frequency. Variation of the linear polar- 
ization position angle within the mean pulse interpreted 
in terms of the rotating vector model (Radhakrishnan & 
Cooke 1969) yields the angle between the observer's line of 
sight and the rotation axis, C — 40°, and the angle between 
the magnetic and rotation axes, a ~ 35° (Manchester & 
Johnston 1995). 

ROSAT observations of PSR J0437-4715 with the Po- 
sition Sensitive Proportional Counter (PSPC) have re- 
vealed (Becker & Triimper 1993; hereafter BT93) that this 
is also a bright (count rate = 0.204 ± 0.006 s'^) soft X- 
ray pulsar with a single broad pulse and a pulsed fraction 
/p = 33 ± 3% in the PSPC energy range 0.1 - 2.4 keV. 
BT93 found that the pulsed fraction varies with photon 
energy E and peaks in the range 0.6 — 1.1 keV, reach- 
ing 53 ± 6%. Although the pulsar spectrum can be fitted 
with a single power law, indicative of a non-thermal ori- 
gin of the soft X-ray radiation, the energy dependence of 
the pulsed fraction in the narrow energy range makes this 
interpretation hardly plausible. BT93 showed also that a 
single blackbody model does not fit the PSPC spectrum 
leaving a residual hard excess above 0.4 keV. They sug- 
gested that the spectrum can consist of two components: 
a power law, representing magnetospheric or nebular emis- 
sion, and a blackbody component of the temperature 
T - 1.7 X 10® K emitted from an area of ~ 0.08 dfgo km^ 
(di80 — d/180 pc); this thermal component was suggested 
by BT93 to be radiated from a hot spot on the neutron 
star (NS) surface. 
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PSR J0437-4715 was also detected with the EUVE 
Deep Survey Instrument (DSI) in the Lcxan band fil- 
ter {E ~ 0.05 - 0.2 keV) by Edelstein, Foster & Bowyer 
(1995) and Halpern, Martin & Marshall (1996; hereafter 
HMM96). According to Edelstein et al. (1995), the DSI 
count rate of the source is 0.0143 ± 0.0008 s~^, whereas 
HMM96 report the count rate of 0.00973 ± 0.00017 s'^ 
obtained with a much longer exposure (496 ks vs. 72 ks). 
Edelstein et al. (1995) ruled out the single power-law spec- 
tral model based on the improbably high hydrogen col- 
umn density, uh = 2.5 x 10^" cm^^, required by this 
model. They claimed that both the ROSAT flux below 
0.4 keV and the EUVE flux could arise from an isother- 
mal blackbody with a temperature ^ 5.7 x 10^ K, an 
emitting area of ~ 3 km^, and an absorbing column of 
nij = 5 X 10^'' cm-2. On the other hand, HMM96 con- 
cluded that the combined analysis of the ROSAT spec- 
trum and EUVE flux is consistent with a single power-law 
spectrum of photon index 7 = 2.2 — 2.5 and intervening 
column density nn = (5 — 8) x 10^^ cm^^. Alternatively, 
the combined data can be interpreted as comprised of two 
components, e. g., a power law and a blackbody compo- 
nent emitted from a hot polar cap of radius 50-600 m and 
temperature (1.0 - 3.3) x 10'' K. HMM96 observed the 
pulsar in the high time resolution mode, which enabled 
them to obtain the light curve and to measure the pulsed 
fraction /p = 27 ± 5% in the EUVE DSI spectral range. 

Thus, although very important data have been col- 
lected and analyzed, the true nature of the soft X-ray ra- 
diation from PSR J0437-4715 remains elusive. The only 
firmly established facts are that the radiation is pulsed 
(pulsed fraction apparently depends on energy), and the 
spectrum cannot be fitted with a single blackbody model. 
The main question is whether the radiation is of a non- 
thermal (magnctosphcric? nebular?) origin or at least a 
fraction of it can be interpreted as thermal (or thermal- 
like) radiation from some heated regions (polar caps?) on 
the NS surface. 

Virtually all the (different) models of radio pulsars 
(e. g., Cheng & Ruderman 1980; Arons 1981; Michel 
1991; Bcskin, Gurcvich & Istomin 1993) predict these 
objects to have polar caps (PCs) around the NS mag- 
netic poles heated up to X-ray temperatures by rela- 
tivistic particles and gamma-quanta impinging onto the 
pole regions from the acceleration zones. A conventional 
assumption about the PC radius is that it is close to 
the radius within which open magnetic field lines origi- 
nate from the NS surface. For PSR J0437-4715, it gives 
iipc = 1.9(i?Ns/10km)3/2 (P/5.75ms)-i/2 km. Expected 
PC temperatures, Tpc 5 x 10^ - 5 x 10^ K, and lumi- 
nosities, Lpc ~ 10-^^ — 10^2 erg s~^, are much less certain, 
being strongly dependent on the specific pulsar model. 
Thus, one cannot firmly predict PC properties because 
of the lack of a well-established pulsar model — rather a 
theoretical model should be chosen based on X-ray obser- 
vations of radio pulsars. 



X-ray observations of other old pulsars (e. g., PSR 

B1929-hl0 - Yancopoulos, Hamilton & Hclfand 1994; PSR 
B0950-t-08 - Manning & Willmore 1994) do show pulsed 
X-ray radiation which, in principle, could be the thermal 
PC radiation. However, the number of photons collected 
has been too small to make firm conclusions, and the oppo- 
site hypothesis, that this radiation is of a magnetospheric 
origin (Ogelman 1995; Becker & Triimpcr 1997 - hereafter 
BT97), cannot be excluded. In principle, one could also 
observe the PC radiation from younger pulsars and use 
these data to discriminate between different PC models. 
Indeed, there are some indications that hard components 
of the soft X-ray spectra of, e. g., PSR B0656-M4 and PSR 
B1055-52 (Greiveldinger et al. 1996) may consist of two 
subcomponents, a power law and a thermal subcompo- 
nent corresponding to emission from PCs of temperatures 
Tpc - 1.5 x 106 K (for PSR B0656-f 14) and - 3.7 x 10^ K 
(for PSR B1055-52). This interpretation, however, is not 
unique because it is difficult to separate the hard compo- 
nent from the soft one which is believed to originate from 
the (cooler) entire NS surface, and even more difficult to 
separate the two subcomponents of the hard spectral tail. 

One cannot also exclude a, priori that the nonther- 
mal radiation from the pulsar magnetosphere or a pulsar- 
powered compact nebula contributes to, or even domi- 
nates, the observed soft X-ray flux of PSR J0437-4715. 
It is natural to assume that the spectrum of this radia- 
tion can be approximated by a power law in the relatively 
short EUVE-ROSAT range. One could also expect the 
magnetospheric (but not nebular) radiation to be pulsed 
with the radio pulsar period; the pulses should be, as a 
rule, narrower than those of the thermal PC radiation, 
and the shape of the light curve should not vary consider- 
ably with photon energy. Thus, if more data confirm the 
conclusion of BT93 that the pulsed fraction depends on E, 
the radiation of PSR J0437-4715 should either be thermal 
or consist of the thermal and nonthermal components. In 
the latter case, since the sources of the thermal and non- 
thermal radiation are expected to be spatially separated, 
it is natural to expect that pulsations of these two com- 
ponents should be phase-shifted due to the difference of 
travel times and aberration. For instance, if the nonther- 
mal component is generated at a distance comparable to 
the light cylinder radius {R\c ~ Pc/2tt ^ 3 x 10^ cm), the 
time delay leads to a phase shift of about 0.2 of the pul- 
sar period. No energy-dependent phase shifts have been 
reported for this source, although it can be explained by 
the poor photon statistics in hard PSPC channels. 

As in the case of thermal PC radiation, the cur- 
rent theoretical models of nonthermal high-energy pul- 
sar emission (e. g., Sturner, Dermer & Michel 1995; Ro- 
mani 1996) are not enough elaborated to predict the 
intensity, spectrum and light curve for a given pulsar. 
There exist empirical formulae (Seward & Wang 1988; 
Ogelman 1995; BT97) which relate the (presumably non- 
thermal) X-ray luminosity to the pulsar parameters, e. g., 
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to the period and magnetic field. For instance, Ogelman 

(1995) pointed out that the observed soft X-ray lumi- 
nosities satisfy, for 7 pulsars, the following equation: 

~ 6.6 X 10^6 (Si2/P2)2-^ erg s'^ oc E^-^^. BT97 fitted 
the luminosities of 26 pulsars observed in the ROSAT 
range with the dependence — O.OOli?. These pulsars 
represent a very wide range of spin-down luminosities 
{W^ -l(y^'^ erg s-i), ages (10-^-7x10^ yr), magnetic field 
strength (10^ — 10^'^ G) and spin periods (1.6 — 530 ms). 
The inferred dependence indicates that for most pulsars 
the bulk of observed radiation is directly connected with 
pulsar mechanisms, i. e., with production and acceleration 
of relativistic particles which carry away the NS rotational 
energy. The fact that X-rays from some of these pulsars 
(e. g., Crab) are certainly nonthermal may allow one to as- 
sume that the radiation detected with the ROSAT from 
most pulsars is of a nonthermal origin. On the other hand, 
a correlation between and E should take place also for 
radiation emitted by the pulsar PCs because their (ther- 
mal) luminosity is also provided by relativistic particles 
generated in the magnetosphere and accelerated towards 
the NS surface, so that the PC luminosity should be a 
fraction of E which goes to heating of the PCs. Moreover, 
both theoretical (see Section 5.1) and observational esti- 
mates of the PC luminosities may happen to be very close 
to the values predicted by the (E) dependence obtained 
by BT97. For instance, the luminosity of the thermal com- 
ponent in the power-law plus blackbody fit for PSR J0437- 
4715, Lboi = 2.2 X 10^0 dfso erg s"! (BT93), is close to the 
predicted value, 4 x 10'^" erg s^^. Another example is PSR 
B1929+10 whose luminosity Lboi ^ 1-2 x 10^" d'^^^ erg s'^ 
inferred from the blackbody fit of both ROSAT and 
ASCA data (Yancopoulos et al. 1994; Wang & Halpern 
1997) perfectly matches the dependence derived by BT97. 
In general case, we may expect that the pulsar X-ray ra- 
diation contains both thermal (PC) and nonthermal com- 
ponents, both growing with E, and the relation between 
the thermal and nonthermal fluxes may be different for 
different objects, depending, in addition to E, on other 
intrinsic pulsar parameters (e. g., pulsar period P, mag- 
netic inclination a), as well as on the rotational inclination 
C (because the radiation beam widths are different for the 
thermal and nonthermal components). Thus, one cannot 
rule out that the PC component dominates in some cases, 
and we explore this possibility for PSR J0437-4715. 

An important evidence on the nature of radiation of 
PSR J0437-4715 could come from deep observations of 
its high-energy tail, at energies above 1 — 2 keV. Although 
this object has been detected by ASCA (two 20 ks obser- 
vations) below 3 keV, the results are still not very con- 
clusive because of poor choice of observing modes which 
did not take into account the presence of a neighbor- 
ing Seyfert galaxy. Nevertheless, preliminary results show 
that a spectrum of ~ 400 pulsar photons extracted from 
the CCD away from the Seyfert galaxy is softer than a 
power law and resembles more a thermal-like spectrum 



(Kawai, Tamura & Saito 1996). This conclusion, however, 
needs verification based on longer pointed observations in 
a mode minimizing contamination from the Seyfert. 

The conclusion of BT93 and HMM96 that the PSPC 
spectrum of PSR J0437 4715 cannot be fitted with a sin- 
gle thermal component is based on the assumption that 
its spectrum coincides with the spectrum of the black- 
body radiation. However, spectra emitted by real bodies, 
including stars, are always different from the Planck spec- 
trum. In particular, if the temperature of a stellar atmo- 
sphere grows inward, and the absorption coefficient de- 
creases with frequency (e. g., fc^ oc for the Kramers 
law), then the spectrum is harder than the blackbody 
spectriim at high frequencies because we see deeper and 
hotter layers. This means that if a NS is covered with a 
fully ionized plasma, whose opacity can be described by 
the Kramers law, its spectriim is substantially harder at 
high energies, hi' > kTes (e. g., Pavlov & Shibanov 1978). 
This general property has been demonstrated in models 
of hydrogen and helium atmospheres of cooling NSs by 
Romani (1987), Pavlov et al. (1995), Rajagopal & Ro- 
mani (1996), Zavlin, Pavlov & Shibanov (1996; hereafter 
ZPS96). One can expect that a NS has a purely hydro- 
gen atmosphere if it experienced accretion of the inter- 
stellar matter — heavy elements of the accreted matter 
sink down rapidly due to the strong NS gravitation (Al- 
cock & Illarionov 1980). Since PSR J0437-4715 is a very 
old object, it is quite plausible that it accreted some mat- 
ter during its long life, and in this case we may expect 
that the excess of the observed flux at E > 0.4 — 0.6 keV 
is due to the fact that PSR J0437-4715, including its PCs, 
is covered by a fully ionized hydrogen atmosphere. 

The PC radiation should be inevitably pulsed unless 
the rotation axis coincides with either the line of sight 
or the magnetic axis. If it were the blackbody radiation, 
the shape of the light curves and the pulsed fraction /p 
would remain the same at all photon energies E. However, 
the atmosphere radiation has a very important feature — 
it is anisotropic, with anisotropy (and consequently light 
curves) depending on energy (Pavlov et al. 1994; Shibanov 
et al. 1995; ZPS96). This dependence is different for dif- 
ferent chemical compositions and surface temperatures. 
E. g., for atmospheres with relatively high temperatures 
~ 10^ K consisting of light elements (hydrogen, helium), 
the anisotropy increases with increasing E in the soft X- 
ray range. This also makes the spectra dependent on the 
rotation phase. The inter-dependence of the spectral and 
angular distributions means that the proper interpreta- 
tion of the observed PC radiation implies fitting of both 
the spectra and the light curves with the same PC model. 

The above-described properties of radiation emitted 
by a PC covered with an atmosphere warrant a new in- 
vestigation of the soft X-ray radiation observed from PSR 
J0437-4715, based on the NS atmosphere models. To per- 
form this investigation, we used the ROSAT and EUVE 
data analyzed previously by BT93 and HMM96 and com- 
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bined with the results of new observations carried out with 

the ROSAT High-Resolution Imager (HRI) and ROSAT 
PSPC (Becker et al. 1997). These observations are briefly 
described in Section 2. 

We show that with allowance for the properties of 
NS atmosphere radiation these joint ROSAT and EUVE 
data can be interpreted as thermal radiation of two PCs, 
without invoking an additional nonthermal component. In 
Section 3 we explore a simplest single-temperature PC 
model which assumes that both PCs have equal radii Rpc 
and uniform temperatures Tpc, whereas the rest part of 
the NS surface has a much lower temperature and is in- 
visible in X-rays. Since the magnetic field of PSR J0437- 
4715 is very low, it cannot affect the radiative properties 
of the PC atmospheres. Therefore, we use the low-field 
atmosphere models developed by ZPS96. To calculate the 
flux as measured by a distant observer, wc integrate the 
specific flux over the visible PC surface with allowance 
for the gravitational redshift and bending of the photon 
trajectories. For the case of small PCs (-Rpc <S; -Rxs) with 
uniform temperatures, a convenient expression for the ob- 
servable flux is given by Zavlin, Pavlov & Shibanov (1995; 
sec their Eq. [A15])). We show that this simple model is 
generally consistent with the observational data if the PCs 
are covered with hydrogen or helium, whereas the iron at- 
mosphere model does not fit the sp{x;tra (cf. Pavlov et 
al. 1996b). For the hydrogen and helium models, we esti- 
mate i?pc, Tpc and the interstellar hydrogen column den- 
sity, Uh- 

The fit with the simplest PC model is still not perfect. 
This is not surprising because it is hard to expect the 
real PCs to be uniformly heated. On the contrary, due to 
higher heat conduction of subphotospheric layers, where 
the energy of accreting relativistic particles is released, the 
heat can propagate along the surface, eventually heating 
surface layers out of the "primary" hot spot. This should 
result in a larger hot region with the temperature decreas- 
ing outward. In fact, this mechanism should be more ef- 
ficient just for low-field pulsars because the strong mag- 
netic fields of ordinary pulsars greatly reduce the trans- 
verse conductivity (e. g., Hernquist 1985). To the best of 
our knowledge, there have been no reliable calculations 
of the temperature distribution around the pulsar mag- 
netic poles. Thus, to include this effect into consideration, 
we explore in Section 4 a simple model: the distribution 
is assumed to be a two-step function ( "core+rim'' ) with 
two temperatures, Tcore and Trim, and two radii, i?core 
and -Rrim- We show that this model is fairly consistent, 
for hydrogen-covered PCs, with all the data available. We 
discuss the results and implications of our interpretation 
in Section 5 and draw conclusions in Section 6. 



2. Observational data and the power-law fit 



2.1. ROSAT PSPC and HRI 

We used two data sets obtained during pointed PSPC ob- 
servations (total 3172 ± 58 counts) and one data set ob- 
tained with HRI (1517 ± 39 counts); the corresponding 
count rates are 0.202 ± 0.004 s'^ and 0.043 ± 0.004 s-\ 
respectively. Details of these observations are described by 
BT93 and Becker et al. (1997). 

Since no source counts were detected in the PSPC 
pulse height channels above channel 200, and the PSPC 
response matrix is not sufficiently known for softest chan- 
nels, we used channels 10-200 {E w 0.1 - 2.0 keV) for 
our spectral and temporal analysis. We manually binned 
the count spectrum into 49 energy bins with the signal- 
to-noise ratio between 4 and 9 and used this binning for 
the spectral fitting with the aid of the MIDAS/EXSAS 
software (Zimmerman et al. 1994). 

We used the PSPC and HRI temporal dependences 
of the source count rate folded with the radio pulsar pe- 
riod by Becker et al. (1997). The PSPC light curves were 
extracted in 0.1-0.6, 0.6-1.1, 1.1-2.0 and 0.1-2.0 keV en- 
ergy ranges (13, 9, 9 and 17 phase bins, respectively). The 
shapes of the light curves from the second data set were 
compared with those of the first set using a test. The 
analysis showed that the light curves from both sets of the 
PSPC data are consistent with each other with a proba- 
bility > 70%, that justifies summation of the light curves 
from the two PSPC observations. 

The combined PSPC light curves show an evidence 
that the pulsed fraction /p depends on energy, in accor- 
dance with the dependence reported by BT93 for the first 
data set. Defining fp in the standard way, as the fraction of 
the total counts lying above the light curve minimum, we 
obtained an increase from /p = 32 ± 4% at 0.1-0.6 keV to 
fp = 44 ±8% at 0.6-1.1 keV. Formally, in the 1.1-2.0 keV 
range the pulsed fraction reaches 54± 11%, but this value 
should be considered with caution because of poor statis- 
tics (~ 252 counts in this range). The pulsed fraction in 
the total energy range is /p = 30 ± 4%. No phase shifts 
between the light curves in different energy ranges were 
found (cf. BT93). 

Although the estimated fp values indicate that the 
light curves have different shapes in different energy 
ranges, they suffer from large uncertainties of the number 
of counts in the phase bins corresponding to the minima of 
the observed light curves, being thereby not very reliable 
estimators. Besides, light curves with equal fp can have 
quite different shapes. Therefore, an independent statis- 
tical test is needed to ascertain and quantify statistically 
the difference of the light curve shapes. The same test 
shows that the probability that the shapes of the light 
curves in different ranges are different is generally high, 
between 77% and 92%. Although the difference might be 
partly due to systematic uncertainties, we believe that the 
dependence of the light curve shapes and pulsed fractions 
on energy is statistically reliable (at least, it cannot be 
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rejected on statistical grounds) and can be considered as 
an additional argument in favor of thermal origin of the 
soft X-ray radiation from PSR J0437-4715. 

The HRI light curve folded with the radio pulsar period 
and binned into 11 phase bins results in the pulsed fraction 
/p = 32 ± 5%. 

2.2. EUVE DSI 

We used the observational data kindly provided by Jules 
Halpern and described in detail by HMM96. The total of 
4163 ±65 source counts were collected at the source count 
rate of 0.00973 ± 0.00017 8"^ This uncertainty (±2%) is 
due to statistics only; a systematic uncertainty may be 
as high as 15-20% (Bowyer et al. 1996). The produced 
light curve with 20 phase bins yields the pulsed fraction 
/p = 27 ± 5%. 
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Fig. 1. 68%, 90% and 99% confidence contours for the 
power-law fit to the PSPC spectrum combined from two data 
sets. Thiclc dot-dashes are the lines of minimum and maximum 
DSI fluxes if only statistical errors are taken into account; thick 
dots correspond to possible ±15% systematic uncertainty of 
the DSI effective axea. 



2.3. Power-law fit 

A power law with the photon index 7 = 2.35±0.20, intrin- 
sic source luminosity = (7.3 ± 0.8) x 10"^° d\^Q erg s^^ 
in 0.1-2.0 keV and = (7.5 ± 2.3) x lO^^ cm'^ pro- 
vides an acceptable (xu = 1-08) fit to the combined PSPC 
spectrum. The parameter values are generally consistent 
with those obtained by BT93 and HMM96: 7 = 2.6 ± 0.2, 
Uh = (1.4±0.5) X 10^0 em-2, = 9.9 x lO^o dj^^ erg s-\ 
and 7 = 2.45 ± 0.25, uh = (9.2 ± 4.8) x lO^^ cm'^, 
Lj, = 8.3 X 10^° dfgo erg s~^, respectively. The best power- 
law fit to the combined PSPC spectrum yields the HRI 



count rate 0.044 s , which is compatible with the mea- 
sured value of 0.043 ± 0.004 

Folding the same model with the EUVE DSI response 
gives the count rate 0.00747 s~^, which is about 20% lower 
than the observed value. Following HMM96; we tested 
whether or not the single power-law fits to the PSPC 
spectrum can be consistent with the DSI flux. We folded 
each of the trial spectra on a grid of the fitting param- 
eters 7 and UH through the effective area curve of DSI 
and obtained a corresponding grid of predicted DSI count 
rates. Figure 1 presents the 68%, 90% and 99% confidence 
contours of the PSPC fits in the j-nH plane, together 
with the bands allowed by estimated systematic uncer- 
tainties of the DSI response. Figure 1 shows that our joint 
PSPC±DSI analysis restricts the hydrogen column den- 
sity to nH = (5.4-7.1) X 10^^ cm~^, at a 90% confidence 
level. This restriction is compatible with that of HMM96, 
riH = (5 — 8) X 10^^ cm~^, obtained with almost thrice 
lower number of the PSPC counts. It should be noted, 
however, that if one adopts the 2% statistical uncertainty 
for the DSI domain, it does not overlap the 68% PSPC 
confidence ellipse, so that the joint analysis of the PSPC 
and DSI data hints that the power-law spectral model may 
be only marginally acceptable. This can be considered as 
one more reason to test thermal models for the soft X-ray 
radiation of PSR J0437-4715. 

3. Single-temperature polar cap model 

In this Section we analyze the above-described data in 
terms of the simplest one-component thermal model. We 
assume that the observed radiation is emitted by two sym- 
metric, uniformly heated PCs around the (dipole) mag- 
netic poles. We consider three chemical compositions for 
the PC atmospheres, and, for the sake of comparison, the 
blackbody model for the PC radiation. The fitting param- 
eters for this model are hh, the PC (effective) temperature 
Tpc, and the ratio Rpc/d. The radius Rpc is defined as the 
distance from the cap edge to the magnetic axis. We fix the 
distance to the pulsar, d = 180 pc (a most probable value, 
according to Sandhu et al. 1997), and present the results 
for Rpc and PC luminosity for this distance. Since the at- 
mosphere spectra and, particularly, light curves depend 
on the viewing angle, the directions of the magnetic and 
rotation axes are to be specified. We adopt, for most of our 
fits, C = 40° for the angle between the rotation axis and 
the line of sight, and a = 35° for the angle between the ro- 
tation and magnetic axes (Manchester & Johnston 1995). 
To specify the gravitational acceleration for calculating 
the atmosphere models, and the NS mass-to-radius ratio 
for account of the gravitational bending of the photon tra- 
jectories, we adopt A/ns = 1-4M0 and i?NS = 10 km for 
the NS mass and radius. All temperatures and radii here 
and below are as measured at the pulsar surface. 




Fig. 2. Single PC model fits (blackbody radiation and iron atmosphere) to the PSPC count rate spectrum, with values of 
best-fit parameters. 



3.1. PSPC spectral fits 

Fitting the combined PSPC spectrum with four models 
for the PC radiation yields the following results. 

Blackbody. — In agreement with what was found by 
BT93 and confirmed later by HMM96 for the first PSPC 
data set, the single blackbody model fails to fit the com- 
bined PSPC spectrum (minimum = 4.78) because the 
blackbody spectrum is too soft (steep) at high energies 
and leaves too large residuals at ii^>0.9keV (Fig. 2). 
Formally, the fit yields uh = (0.0 ± 1.1) x lO^'' cm'^, 
logTpc = 6.28 ± 0.03 and i?pc = (0.17 ± 0.02) diso km. 
This gives the bolometric luminosity of two PCs Lboi = 
(1.36 ± 0.25) X 10^° dlso erg 8"^ 

Iron atmosphere. — Because of numerous spectral lines 
quasi-randomly distrubuted at ~ 0.3 — 3.0 keV, the 
smoothed spectra emitted by iron atmospheres in the soft 
X-ray range are not much different from blackbody spec- 
tra at the same temperatures (Rajagopal & Romani 1996; 
ZPS96). The spectral lines of iron atmosphere spectra can- 
not be detected by ROSAT due to its poor energy reso- 
lution. The spectral fit with the iron atmosphere model 
(right panel of Fig. 2) is even worse than with the black- 
body: minimum xt = 10.46. The lower quality of this fit 
is caused by the fact (ZPS96) that the iron atmosphere 
spectra are softer the blackbody spectra in the soft X-ray 
range at temperatures ^ 10^ K. 

Hydrogen atmosphere. — In comparison with the two 
cases above, the PC model with hydrogen composition 
gives a much better fit, with minimum xt = 1-68 (Fig. 3). 



Since the hydrogen atmosphere spectra, as well as any 
other thermal spectra, have an intrinsic decrease towards 
low energies, they fit the observed spectrum at much lower 
interstellar hydrogen density, Hh = (0.0 ± 1.7) x 10^^ 
cm~^. The other parameters are logTpc = 5.94 ± 0.03 
and Rpc = (0.77 ± 0.13) diso km. Since the high-energy 
tail of the radiation emitted from hydrogen atmospheres 
is harder than the blackbody tail, the resulting PC tem- 
perature is 2.2 times lower than the temperature of the 
blackbody fit. To provide the same fltix, Rpc becomes 
4.5 times greater approaching the conventional estimate 
of ^ 1.9 km. The corresponding bolometric hmiinosity of 
two PCs is iboi = (1-22 ± 0.22) x 10^° dl^g erg s-\ com- 
parable with that for the blackbody model. We checked 
that the fitting parameters are very stable with respect 
to the choice of energy channels (binning) of the PSPC 
count rate spectrum — only the minimum value of xt 
varies slightly. For example, the fit in 20-140 channels 
(0.2 — 1.4 keV range) with 19 bins yields the same val- 
ues of Uh, Tpc and Rpc, with minimum xt = 1-16. 

Helium atmosphere. — For helium- covered PCs (right 
panel of Fig. 3), the quality of the fit and the values of the 
model parameters are similar to those for the hydrogen 
composition. A smaller value of the minimum Xv = 1-37 
is due to the fact that at T ^ 10^ K the spectrum of 
the helium atmosphere is slightly harder than that of the 
hydrogen atmosphere (ZPS96). Generally, the fits with hy- 
drogen and helium models are almost indiscernible. 




Fig. 3. The same as in Fig. 2 for the hydrogen and heUum atmospheres. 



Although the hydrogen and hehum PC models fit 
the PSPC spectrum much better than the blackbody 
model, there still is an excess of the observed radiation 
at E > 1.5 keV (see Fig. 3). The excess can be partly due 
to systematic errors of the high-energy PSPC response 
which may be as high as 20%. Irrespective of that, the ex- 
cess can be naturally explained by simplicity of the single- 
temperature PC model — if there is a smaller area of a 
greater temperature within the PC, it can easily provide 
the excess high-energy quanta, the possibility we explore 
in Section 4. The single-temperature approximation, how- 
ever, enables us to show that the spectrum can be fitted 
with a purely thermal model provided that the PCs are 
covered by a hydrogen or helium atmosphere, whereas the 
hypothesis about the iron surface of the PCs can be firmly 
rejected. (A similar conc;Iusion has been drawn by Ra- 
jagopal & Romani (1996) who, however, obtained substan- 
tially different PC parameters because they fitted the ob- 
served spectrum without allowance for the limb-darkening 
and light bending.) We do not expect that other heavy- 
element compositions of the NS surface could fit the ob- 
served spectrum — they should also lead to emitted spec- 
tra much softer that the spectra of hydrogen or helium 
atmospheres. We consider the hydrogen atmosphere more 
natural that the helium one (e. g., if the atmosphere was 
formed due to accretion of the hydrogen-rich interstellar 
matter or matter from the binary companion) , and we ex- 
plore consistency of this hypothesis with the HRI and DSI 
results, as well as with the light curves observed with the 
three instruments. 



3.2. Analysis of the ROSAT HRI and EUVE DSI count 
rates 

Figure 4 shows the confidence contours for the parameters 
of the hydrogen PC fit to the PSPC spectrum. We checked 
that any point from the 99% PSPC confidence region pro- 
duces the model HRI count rate within the observed range. 
The best parameters of the PSPC spectral fit give a too 
high count rate for the EUVE DSI, 0.01136 s-\ which 
does not even get into the ±15% systematic uncertainty 
range around the observed value. The reason is that the 
formal best fit yields uh = (no interstellar absorption), 
which clearly is not a realistic value. However, nn is deter- 
mined with some uncertainties, and even a small value of 
uh can noticeably decrease the DSI flux and make the DSI 
and PSPC count rates consistent with each other. Similar 
to the case of the power- law model (Section 2), we com- 
puted the bands in the parameter planes allowed by the 
observed DSI fiux. The thick dot-dash lines in Fig. 4 corre- 
spond to the DSI count rates of 0.00973 ± 0.00017 cnts/s. 
The shaded regions show the domains of the model pa- 
rameters allowed by all the three instruments at the 68% 
confidence level. If we adopt the ±15% uncertainty for the 
observed DSI count rate (thick dots), the allowed regions 
become much broader. 

Thus, there exist a domain of the model parame- 
ters compatible with all of the detected count rates and 
with the PSPC spectra: uh = (0.5 - 0.9) x lO^'' cm"^ 
logTpc = 5.92 - 5.95 and Rp^ = (0.75 - 0.85) diso km. 
The corresponding area and bolometric luminosity of two 
PCs are A = (3.5 - 4.5) dfgo km^ and iboi = (1-2 - 
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Fig. 4. Projections of the 68, 90 and 99% confidence domains onto two parameter plains for the hydrogen PC fit to the PSPC 

spectrum. Thin dashes correspond to constant values (numbers near the lines) of third parameter. Filled squares show the 
best-PSPC-fit parameters (cf. Fig. 3). Thick dot-dashes and dots are the lines of constant DSI fluxes (cf. Fig. 1). The shaded 
regions give the model parameters allowed by the count rates of all the three instruments. Crosses correspond to the "best PC 
model". 



1.3) X 10^° dfgQ erg s~^. Similar to the case of the power- 
law fit, the main role of the EUVE data is that these 
data enable us to restrict the hh value; the difference is 
that for the thermal PC model the range of jointly ac- 
ceptable nn values lies above the best-PSPC-fit value, 
and they are an order of magnitude lower than for the 
joint power-law fit. This implies a higher ionized frac- 
tion, and a lower mean number density of the interstellar 
hydrogen in the direction to PSR J0437-4715 (see Dis- 
cussion). The center of the allowed three-dimensional pa- 
rameter domain is given, approximately, by the follow- 
ing values: hh = 0.7 x 10^^ cm~^, logTpc = 5.93 and 
Rpc = 0.82 di80 km, which can be considered as a best 
joint PSPS-hHRI+DSI fit (marked with crosses in Fig. 4). 
For the sake of brevity, we call this set of parameters "i/ie 
best PC modeP for the rest of this Section. 



3.3. Light curve analysis 

Now we should test whether the PC model obtained from 
the analysis of the PSPC spectrum and HRI and DSI 
count rates fits the light curves detected with these in- 
struments. To do that, we computed the spectra for 80 
rotation phases, < < 1, for the best PC model and 
convolved the spectra at each phase with the PSPC, HRI 
and DSI responses. Then, we binned these model count 
rates into the the phase bins adopted for the observed 
light curves to compare the light curves with the aid of 
the x^-test. The value depends on how we co-align 
the phases for the model and the observed light curves. 
Since the reference phases of the observed light curves 
were chosen arbitrarily, we may vary the phase shift be- 
tween the model and observed light curves to minimize x^- 
(Of course, the shifts in different energy ranges should not 
differ from each other by more than a fraction of the bin 
width; ideally, the difference should be comparable with, 
or less than, the instrument time resolution). 



Figure 5 presents the adjusted model and observed 
PSPC light curves for the four energy ranges. We de- 
fine the zero rotation phase as the phase when the an- 
gle 6 between the magnetic axis and the line of sight 
equals its minimum value, ^min = Q — a = this phase 
corresponds to the maximum of the model light curve, 
whereas the minimum lies at (j) = 0.5, corresponding to 
~ Smax ~ C + ct ~ 75°. The difference of the phase shifts 
applied to minimize x^ does not exceed 1/80 of the rota- 
tional period, which corresponds to 0.07 ms, about one 
half of the the ROSAT time resolution. As is seen in 
Fig. 5, the best PC model yields the light curves in ex- 
cellent agreement with those observed in the 0.1 — 0.6 and 
0.6 — 1.1 keV ranges not only the pulsed fractions arc 
within the statistical uncertainties, but also the xt values 
are small enough. For the 1.1—2.0 keV range, the model /p 
is within the observed uncertainty, but the fit is formally 
unacceptable {xt = 3.81). This is not surprising in view 
of the aforementioned deficit of PC counts with respect to 
observed in this energy range. 

Figure 6 shows the observed ROSAT HRI and EUVE 
DSI light curves It is worth noting that the difference be- 
tween the phase shift applied to adjust the observed PSPC 
and HRI light curves to the model ones is as small as 1/80. 
The coincidence of the model and observed light curves is 
fairly good for both HRI and DSI. 

The results above were obtained for ( = 40°, a = 35°, 
and the standard NS mass and radius. The angles, how- 
ever, are not very certain as this pulsar shows a com- 
plicated variation of linear polarization angle across the 
eight-component mean radio pulse (Manchester & John- 
ston 1995). In principle, the interpretation of the soft X- 
ray data adopted in the present paper can reduce this un- 
certainty — one can fit the data with the PC model consid- 
ering a and C as fitting parameters. These angles, as well 
as the (unknown) ratio Mns/-Rns, affect insignificantly 
the spectral fit parameters {uh, Tpc and Rpc). However, 
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Fig. 6. The same as in Fig. 5 for the ROSAT HRI and EUVE 
DSL 
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Fig. 7. Two examples of hydrogen PC light curves: for another 

sot of the angles {( = 24° and a — 20°; solid) and for the 
best PC model without allowance for the gravitational bending 
(dashes) . 



Fig. 5. Observed and model PSPC light curves in four energy 
ranges. The model curves correspond to the "best PC model" 
(model parameters in the upper panel) . 



the shape of the light curves and /p are very sensitive to 
the choice of a, ( and Mns/^^ns- For example, an alterna- 
tive set of the angles, ( = 24° and a. = 20°, obtained with 
another treatment of polarization data (Gil & Krawczyk 
1997), results, for Mns/-Rns = 0.14 Mo/km, in a too low 
pulsed fraction, /p = 15.3% in the whole PSPC energy 
range, so that the PC model is in clear disagreement with 
observations (xt = 3.52) — see Fig. 7. To provide an ac- 
ceptable fit at this set of angles, the pulsar mass is to be 
very low, Mns < 0-5Mq at any i?NS allowed by the equa- 
tions of state of NS interiors (Pavlov et al. 1997). One 
more example shown in Fig. 7 demonstrates the impor- 
tance of the gravitational bending of photon trajectories: 
if this effect is neglected, the best PC model yields a light 
curve with strong pulsations, fp = 76.8%, resulting in a 



great deviation, xt = 29.1. The mass-to-radius ratio can 
be also considered as a fitting parameter to provide better 
coincidence of the light curves and put constraints on the 
equation of state of the NS interiors. 

4. Two-temperature polar cap model 

Motivated by the results of Section 3, we examine a more 
complicated PC model which allows, in the first approxi- 
mation, for a non-uniform temperature distribution along 
the PC surface. We adopt a two-step approximation, 
"core-|-rim" , with Tcore > ^rim and Rcoie < -Rrim- The 
definitions of i?corc and i?iiin are analogous to Rpc in Sec- 
tion 3. This model contains, together with uh, five fitting 
parameters. It is natural to expect that Trim < Tp^ < T^ore 



where Tpc and i?pc arc the pa- 



and i?c 

rameters of the one-component PC model, so that the rim 
and the core should be responsible for the soft and hard 
spectral components of this two-component model. We re- 
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Energy (keV) 



Fig. 8. Two-component ( "core-|-rim" ) fit to the PSPC spectrum with five free parameters (best values in the upper- left panel). 
Photon spectral fiuxes are shown in the right panel for the five-parameter fit (solid curves) and for the fit with fixed values of 
Tcore and -Rcore (dashes). 



tain the same assumptions on the pulsar angles, mass, ra- 
dius and distance as earlier. 



4-1- Iron atmosphere 

As one could expect for a two-component model, the 

iron atmosphere model yields an acceptable fit of the 
PSPC spectrum (minimum xl = 0.72, i' = 44), with the 
best-fit parameters uh = 7.92 x 10^^ cm~^, log Tcore = 
6.57, log Trim = 5.67, i?coro = 0.04 diso km and i?rim = 
8.36 diso km. However, the value of i?rim is very close to 
the adopted pulsar radius i?NS = 10 km. This means that 
the soft component is emitted from almost the whole NS 
surface with the uniform temperature Trim. By virtue of 
this, the rim cannot produce the pulsed fraction fp ~ 30% 
observed in the 0.1-0.6 keV range (the model fp does not 
exceed a few percent). Since the light curves cannot be fit- 
ted with the same model as the spectra, we conclude that 
the hypothesis about the iron NS surface is incompatible 
with the thermal model for the soft X-ray radiation of PSR 
J0437-4715. Thus, we exclude the iron atmosphere model 
from further consideration, as well as the two-component 
blackbody model which also leads to i?rim — -Rns- 



4-2. Hydrogen atmosphere 

The model count rate spectra in Fig. 8 correspond to the 

best spectral fit with the following parameters: nn = 

1.55 X 1019 Cm-2, log Tcore = 6.21, log Trim = 5.63, i?core = 

0.17(^180 km and i?rim = 2.36 digo km (minimum xt = 
0.68, v = 44). The corresponding bolometric luminosity of 
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Fig. 9. Confidence contours for the parameters of hard (core) 
component in the two-component PSPC spectral fit. 



(51% 
is close 



two rims and cores, Tboi = 1-37 x lO^^d^gQ erg s~ 
and 49% from the cores and rims, respectively), 
to that of the single-temperature PC model. 

However, in spite of quite realistic values of these pa- 
rameters, they are very uncertain, which follows from large 
formal errors provided by the fitting code. Although these 
errors, calculated from the derivatives of the fit statis- 
tics, are not reliable, they indicate that the parameters are 
strongly correlated, and the best fit gives only qualitative 
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Fig. 10. Projections of 68, 90 and 99% confidence domains for the three-parameter PSPC spectral fit of the soft (rim) component 

(thin solid and dot lines), with fixed best core parameters (cf. Fig. 9). Thick dashed lines show restriction (68% confidence level) 
obtained from the light curves in the whole PSPC energy range. Other notations are the same as in Fig. 4. 
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Fig. 11. The same as in Fig. 5 for the "core+rim" 
two-component fit. The dashes show the the contribution to 
the o. 1-0.6 keV light curve from the rim component 



estimates of the parameter values. In principle, the param- 
eter ranges can be reduced with the use of the ROSAT 
HRI and EUVE DSI data, but it would require an ex- 
tremely time-consuming analysis of five-dimensional con- 
fidence volumes. A practical way to solve the problem is 
prompted by the shapes of the best-fit component spectra 
— we see that the hot core component gives the main con- 
tribution at higher energies, above 0.4-0.6 keV, whereas 
the rims are responsible for the fiux at lower energies. 
Therefore, we can first fit the hard channels separately 
with a single-temperature cap model (notice that nn of 
expected magnitudes does not affect spectra at these en- 
ergies), and then fix some set(s) of Tcore and i?core and fit 
the rim component with the three fitting parameters, nH, 

Trim and i?rim- 

The first step with 21 energy bins in the 53-200 
channels gives the best-fit parameters of the hot core, 
logTcore = 6.16 ± 0.07 and i?core = (0.23 ± 0.09) diso km 
ixi = 0.92, v = 19), in accordance with the two- 
component results. Figure 9 demonstrates the confidence 
contours for this fit. 

At the second step, we fixed the best pair of Tcore and 
i?core obtained at the first step and fitted the total PSPC 
count rate spectrum with the "core-|-rim" model (the best- 
fit parameters are given in the right panel of Fig. 8). We 
see that the separate fitting of the hard and soft (core 
and rim) components yields spectra fairly close to those 
obtained with the five-parameter fit. 

The projections of the three-dimensional confidence re- 
gions for the parameters of the soft component are shown 
in Fig. 10. The regions are still too broad, in spite of the 
reduced number of the fitting parameters, because the pa- 
rameters are strongly correlated (e. g, Tij{ cx In Rriui at a 
given Trim) due to narrowness of the soft-energy range. In 
particular, the 99% confidence surface allows the outer rim 
radius to be comparable with i?NS at lower Trim and higher 
uh- However, large values of .Rrim are not consistent with 
the strong time modulation of the observed flux at low en- 
ergies, and upper limits for Rrim and Hh (and a lower limit 
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on Trim) are to be obtained from the light curve analysis. 
Since the values of i?rim involved in the light curve fit- 
ting are relatively large, we have to integrate over the rim 
surface making use of the exact equations for the gravita- 
tional bending (Zavlin et al. 1995), instead of previously 
used approximate approach for i?rim Rns- Fitting the 
PSPC light curve in the whole energy range with the two- 
component model restricts the radius: i?rim < 6.2, 6.7 and 
7.1 km; temperature: log Trim > 5.45, 5.42 and 5.39; and 
hydrogen column density: hh/{W^^ cm~^) < 4.4, 4.8 and 
5.2; at 68%, 90% and 99% confidence levels, respectively. 
The projections of the 68% confidence region are shown 
in Fig. 10. 

Further constraints on the fitting parameters follow 
from the comparison with the EUVE DSI count rate (HRI 
count rate is consistent with the PSPC fits and does not 
provide additional restrictions of the fitting parameters). 
The DSI count rate corresponding to the best PSPC fit 
(with the hard and soft components fitted separately) is 
0.01196 s"-*^, i. e., 23% higher than the observed value. We 
used the same approach as for the single-temperature PC 
model and found the three-dimentional volumes in the pa- 
rameter space compatible with the measured DSI count 
rates with the 2% (statistical) and 15% (systematic) un- 
certainties (thick dash-dot and dot curves in Fig. 10). The 
shaded areas in Fig. 10 delimit projections of the parame- 
ter volume compatible with the data collected with all the 
three instruments at the 68% confidence level (neglecting 
poorly-known systematic uncertainties) . We chose a point 
with uh = 3.6 X 10^^ cm~'^, log Trim = 5.53 and i?rim = 
4.0di8o km {^^best core+rim modei'). The bolometric lu- 
minosity for this model is Lboi = 1-57 x lO^^d^gQ erg s"-'^. 

Following the approach used for the one-component 
PC model, we plotted the two-component light curves cor- 
responding to this model (Figs. 11 and 12). The compari- 
son shows a satisfactory agreement with the observed light 
curves. We see that the pulsations at lower energies are 
stronger for the "core+rim" model than for the single- 
temperature PC model (/p = 33.7% vs. 29.0% for the 
0.1-0.6 keV PSPC range, and 25.7% vs. 22.9% for DSI). 
This is due to the fact that the rim, which mainly con- 
tributes a,t E < 0.5 keV (Fig. 8), has a lower temperature, 
Trim < Tpc, resulting in higher anisotropy of the emergent 
radiation (ZPS96). The rim emits 74%' of all the PSPC 
counts with /p = 37.1% at = 0.1 - 0.6 keV, and 87% of 
counts with /p = 31.6% in the DS range. The soft com- 
ponent gives less than a half of the total HRI counts, and 
we obtain almost the same light curve fits for the single- 
temperature PC and "core-|-rim" models (cf. Figs. 6 and 
12). In the whole PSPC range there is a bigger share of 
the total radiation with /p = 36.4% coming from the rim, 
that leads to higher pulsations in the summed light curve 
than for the simplest PC model (/p = 35.0% vs. 32.7%, 
respectively) . 

Finally, we can conclude that the two-component hy- 
drogen "core+rim" model satisfactorily matches to the 
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Fig. 12. The same as in Fig. 11 for the ROSAT HRI and 
EUVE DSI, with separate contribution from the rim compo- 
nent (dashes). 

spectral data and the light curves provided by all the three 
instruments. 

4-3. Helium atmosphere 

As expected, this model yields virtually the same re- 
sults as the hydrogen atmosphere model. For instance, 
the best five-parameter fit gives uh = 1.98 x 10^^ cm~^, 
logTcorc = 6.14, logTrim = 5.57, i?corc = 0.26 di8o km, 
and i?rim = 2.94 diso km, with xt = 0.70. Since, in ad- 
dition, the hydrogen composition looks more natural, the 
helium atmosphere model does not require a detailed sep- 
arate analysis. 

5. Discussion 

5.1. Comparison with predictions of radio pulsar models 

It is interesting to compare the inferred PC properties 

with those predicted by different radio pulsar models. The 
slot-gap model (Arons & Scharlemann 1979; Arons 1981) 
predicts the following energy fiux of relativistic particles 
incident on the stellar surface (Arons 1981) 

/10km\''/' / sina \ 

[-R^) [^) '~ ' 

where ^26 is the magnetic moment in units of 10^^ G cm'\ 
fp is the the ratio of the dipole radius of curvature to 
the actual radius of curvature, 1 $ $ 7 for the parame- 
ters of PSR J0437-4715. The model is applicable for P < 
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Z.lRiQ^^"^ fl^^^"^ ms; this corresponds to fp > 2.85i??o^ for 
P = 5.75 ms, and wc adopt an intermediate value fp = A 
for estimates. If most of this energy is reradiated, the ef- 
fective temperature and luminosity are 
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where fa is the ratio of the polar cap area to its canonical 
value, 11.4i?fQ(5.75 ms/P) km^ [fa ~ 0.5 for the original 
slot-gap model which implies that the particle accelera- 
tion zone is associated with "favorably curved" open field 
lines) . Wc sec that the predicted temperature is 20 — 30% 
lower than the temperature Tpc inferred from the single- 
temperature PC model, but it is higher than Trim- The 
predicted luminosity is in excellent agreement, for /p = 4 
and fa = 0.5, with the observed value. Since the observed 
luminosity is almost independent of the PC temperature 
distribution (2L = (1.0 - 1.5) x 10^° erg s"! for both the 
single-temperature and two-temperature PC models), it 
is the most suitable parameter to compare with, and we 
should conclude that the slot-gap model agrees very well 
with the suggested interpretation of the soft X-ray radia- 
tion from PSR J0437-4715. 

According to the outer-gap model (e. g., Cheng, Ho 
& Ruderman 1986), PSR J0437-4715 is near the pulsar 
"death line", Pdeath = 5.8(5/3.4 x 10^ Gf^^ ms, that 
implies high efficiency of the gamma-ray radiation pro- 
duced by primary e^ accelerated in vacuum gaps in the 
outer magnetosphcre of the pulsar. The total (maximum) 
flux of relativistic e+ (or e~) impinging on the PC from 
the starward end of the acceleration zone can be estimated 



as 



N 



'2ec 



(4) 



where rmin is the minimum distance of the accelera- 
tor from the NS center. Adopting rmin — (l/3)7'ic — 
90(P/5.75 ms) km, the first term in the square brackets 
can be neglected if 7i » 3 x 10^ (the latter value is much 
lower than 7max ^ 10^ — 10* limited by the curvature 
radiation losses). In this case, 7/~2.8xl0®, and the ex- 
pected PC luminosity, icap ~ 3.2 x 10^^ erg s~^, exceeds 
the observed value by a factor of 40-60. The model PC 
luminosity would agree with the observed one only if we 
assume 7^ ~ 7/ ~ 6 x 10^ ~ 10~^7max- Thus, we should 
conclude that either the energies at the starward ac- 
celerator end are much lower than those expected from 
the curvature radiation losses, or there exists a mecha- 
nism reducing the particle flux N which may reach the 
NS surface. In particular, the estimate (5) assumes the 
relativistic particles impinge onto the canonical PC area 
A = (irnR^/c) ~ 11 km^. If the actual area is 40 - 60 
times smaller (a "primary" cap radius is 250 — 300 m, 
compatible with our Pcore), then, at the same (Goldreich- 
Julian) current density, the estimated luminosity coin- 
cides with that observed. An evidence of the reduced effi- 
ciency of the particle accelerator is provided by the upper 
limit on the luminosity of gamma-rays above 100 MeV, 

< 1.7 X 10^2 ^ o.04£; (Fierro et al. 1995). Since the 
same particle accelerator is responsible for PC heating, 
HMM96 suggest that an upper limit for the luminosity of 
one PC can be evaluated as 0.04 times the above-estimated 
maximum value, Leap < 1-3 x 10^" erg s""'^, very comfort- 
able with the observed range of (0.6 — 0.8) x 10^" erg s~^. 

In the model of Beskin, Gurevich and Istomin (1993), 
the energy of relativistic particles impinging onto the PCs 
per unit time is 



surf 



l + K 



E, 



(6) 



where Q = 0.52(P/5.75 ms)"/iO(P/1.99 x 10^^0)^2/5 
K is a. multiplication coefficient, the number of particles 
which are knocked out from the NS surface by an impact- 
ing relativistic particle. The temperature of the polar cap 
is estimated as 



3.8 X 10*^ 



P 



1/20 



(1 -I- xy/^ \ 5.75 ms 1.99- 



20 



K 



(7) 



and the luminosity is Leap = NjfmeC^, where 7/ is the 
Lorentz factor of particles when they reach the NS surface. 
While traveling to the surface, the particles loose a frac- 
tion of their initial energy 7imeC^ via curvature radiation, 
so that (cf. Halpern & Ruderman 1993; HMM96) 



7/ 



1 



3.65 X 106 



1 

3 



In 



RnsJ 



5.65 ms\ / rmiiiN 



-1/3 



-1-1/3 



(5) 



This temperature is consistent with the observed values, 
1 X 10^, only for large multiplication coefficient, K ~ 
200, much greater than the maximum values, JsTmax — 1 
and ~ 20 calculated by Bogovalov & Kotov (1989) for 
7/ = 2 X 10^ and 10^, plausible Lorentz factors for PSR 
J0437-4715. Moreover, the same calculations show that K 
sharply decreases with decreasing B from ~ 10^^ G and 
becomes <C 1 at P < 2 x 10^^ G. Thus, we have to conclude 
that either this model strongly overestimates the PC tem- 
perature and luminosity, or the multiplication coeflacient 
is much greater than calculated. 
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5.2. Implications for the magnetic field geometry and the 
NS mass and radius 

The fits of the spectra and of the hght curves presented 
in this paper were obtained for fixed values of the an- 
gles ( = 40° and a = 35° inferred by Manchester & 
Johnston (1995) from the phase dependence of the po- 
sition angle of the radio polarization. Since these angles 
were determined approximately, and it is not trivial to 
estimate their uncertainties from the radio data, we can- 
not exclude that their true values may be different. We 
also adopted canonical values for the NS mass and ra- 
dius, Mns = 1.4A'/q and i?NS = 10 km, whereas the 
true mass and radius may differ from these values. We 
have demonstrated (see Fig. 7) that the shape of the light 
curves is very sensitive to the values of (, a and M-^g/ R^g. 
This opens a new opportunity to constrain these quan- 
tities using them as free parameters for the light curve 
fitting. This analysis is presented elsewhere (Pavlov et 
al. 1997). Here we only mention that for any allowed val- 
ues of a and (, the mass-to-radius ratio is constrained as 
Mns < 1.7M0(i?Ns/lO km), whereas for a = 35°, C = 40° 
we obtained 1.4 < (MNs/MQ)/(i?Ns/10 km) < 1.6. The 
latter inequality means that if the NS mass equals 1.4M0, 
its radius is 8.8 < R^^s < 10.0 km. 

5.3. Implications for the interstellar hydrogen 

We have shown that interpretation of the soft X-ray data 
in terms of thermal-like radiation from PCs yields system- 
atically lower values for the intervening hydrogen column 
density than the power-law fits of the spectrum. For in- 
stance, BT93 found uh = (1.4 ± 0.5) x 10^° cm'^ for 
the one-component power-law fit, whereas we obtained 
Uh < 2 X 10^^ cm~^ and uh = (1 — 5) x 10^'^ cm~^ for the 
single- and two-temperature PC models, respectively. It 
is useful to compare these values with those inferred from 
observations of other objects close towards PSR J0437 
4715 (/ = 253°, b = -42°; d = 180 pc). The pulsar is in 
the southern third galactic quadrant (180° < I < 270°, 
b < 0°), apparently within the famous low HI column re- 
gion that extends out to 200-300 pc at 210° < I < 265° (at 
least for b > -30° — see, e. g., Paresce 1984; Welsh 1991). 
For instance, two stars of close longitudes and distances, 
eCMa {I = 240°, b = -11°, d = 187 pc) and /? CMa 
{I = 226°, 6 = -14°, = 220 pc), show the neutral hydro- 
gen densities as low as log n/// = 18.0— 18. 5 and 18.2-18.4, 
respectively (Fruscione et al. 1994; Welsh 1991). The den- 
sities in directions of two other stars within the low density 
region, RE 0503-289 (/ = 231°, b = -35°, d = 90 pc) and 
RE 0457-281 {I = 229°, b = -36°, d = 90 pc), with lati- 
tudes closer to that of the pulsar (farther from the Galac- 
tic plane), albeit with twice smaller distances, are even 
lower: lognni = 17.3 - 17.8 and 17.1-17.9. A white dwarf 
WD 0320-540 {I = 267°, b = -52°, d = 103 pc), which is 
likely near the boundary of the low-density region, shows 



lognjj/ = 18.9 — 19.2. Based on these values and observa- 
tions of other (scarce) stars in this region of the sky (Frus- 
cione et al. 1994, Welsh 1991; and references therein), we 
expect that plausible neutral hydrogen column for PSR 
J0437-4715 should not exceed uhi ~ (1 - 3) x 10^^ cm^^^ 
Interstellar absorption of X-rays enters the spectral 
fits through the attenuation function, exp(— njjo-), where 
(7 is the effective cross section per hydrogen atom calcu- 
lated under the assumption that the interstellar hydro- 
gen is not ionized (e. g., Morrison & McCammon 1983). 
The effective hydrogen column density nn coincides, nat- 
urally, with nni for the nonionized ISM, but for the par- 
tially ionized hydrogen we have Hhi < nu < , where 
j^tot ^ _^ jg ^Yie total (neutral plus ionized) hy- 
drogen column density. With the electron column density 
He = 0.81 x 10^^ cm~^, known from the pulsar dispersion 
measure, we can restrict the mean ionization degree of hy- 
drogen, ^ = Ue/nJ^*, as ne/{ne + nu) < i < ne/nu- If we 
adopt Uh = (1—3) x 10^^ cm~^ as a plausible estimate, the 
ionization degree 0.2 < ^ < 0.8 is noticeably greater, and 
the (total) mean number density in the direction towards 
PSR J0437-4715, n^g^d = (0.02-0.07)d^8o cm'^, is lower 
than the estimates, ^ ~ 0.04 - 0.1 and n^^^d ~ (0.16 - 
0.34)d]~gQ cm~'^, obtained from the power-law fit. The in- 
ferred high ionization is in line with ^ « 0.9 found by Gry, 
York & Vidal-Madjar (1985) in the direction to /? CMa. 
It is worth noting that the mean number density inferred 
from our PC model fits is fairly close to the ISM density 
of the ambient medium around the pulsar estimated by 
equating the ram pressure of the interstellar medium to 
the pulsar wind pressure at the apex of the pulsar wind's 
bow shock (see, e. g., HMM96): n*§^ = E / (Anrlv'^cmH) ^ 
0.08 cm~'^, where the distance from the pulsar to the apex 
of the bow shock, ry^ = 2.4 x 10^^ cm, and the velocity of 
the pulsar, Vp = 118 km s~^, are evaluated for d = 180 pc. 

5.4- X-rays from other millisecond pulsars: thermal or 
nonthermal? 

The fact that all available X-ray and EUV data are consis- 
tent with our interpretation still docs not prove that this 
interpretation is unique — it often happens that the same 
data can be equally well fitted with quite different mod- 
els. Therefore, it is important to verify that the hypothesis 
does not contradict to observations of similar objects. Ac- 
cording to BT97, soft X-ray radiation has been observed 
from seven millisecond pulsars. At least one of them, PSR 
B1821-24 (P = 3.05 ms, E = 2.2 x 10^^ erg s'^), shows 
a strong evidence that its radiation is of a nonthermal 
(magnetospheric) origin — its spectrum, observed with 
ASCA up to 10 keV, can be fitted with a power law, and 
its X-ray pulses are sharp (Saito et al. 1997). Thus, if we 
assume that PCs is a phenomenon common for all mil- 
lisecond pulsars, we must explain why this phenomenon 
is not observed for PSR B1821-24. This can be done with 
the use of the aforementioned pulsar models. If the correct 
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estimate obtained for the PC luminosity from the slot-gap 

model is not a chance coincidence, we may expect that the 
luminosity is scaled as Leap oc pi/2p-23/8^ g^j^^ ^^le bolo- 
metric luminosity from two PCs of PSR B1821-24 can be 
estimated as _Lboi ^ (6 — 9) x lO^'^ erg s~^. This estimate 
is much lower than the luminosity = 1.7 x 10^"^ erg 
inferred by Saito et al. (1997). If we adopt the outer- 
gap model and assume that the accelerator efficiency is 
reduced by the same factor for all pulsars, then the lu- 
minosity is scaled as Leap oc pi/2p-ii/6^ which yields 
Lho\ (3 - 5) X lO^i erg s^^ for PSR B1821-24, even less 
than for the slot-gap model. Moreover, even the maximum 
possible PC luminosity predicted by the outer-gap model 
remains lower than observed. Thus, it is not surprising 
that the PC radiation is not seen from this pulsar, being 
overwhelmed by a nonthermal soft X-ray radiation. Saito 
et al. (1997) speculate that the high luminosity of the non- 
thermal radiation may be associated with large values of 
the magnetic field at the light cylinder, Bic ~ 1.5 x 10^ G, 
close to that of the Crab pulsar and much higher than 
Ae ~ 1.7 X 10'' G for PSR J0437-4715. 

The nature of X-ray radiation from other five millisec- 
ond pulsars is much less certain because of small num- 
bers of counts detected. Their luminosities in the ROSAT 
range, inferred under the assumption of a power-law spec- 
trum with the photon index 7 = 2, obey the same depen- 
dence, Lx — O.OOIL^, as ordinary pulsars (BT97), and they 
are greater than the PC luminosities predicted for these 
objects by the slot-gap model (scaled to PSR J0437-4715). 
However, the ratios of the observed to predicted luminosi- 
ties, ranging from ~ 1.6 for PSR J1012-h5307 to ~ 10 for 
PSR J0751-hl907, are smaller than ~ 25 for PSR B1821- 
24. Since the power-law model always yields a higher lumi- 
nosity than the PC model (because of higher hr needed to 
suppress the intrinsically high low-energy flux), one may 
expect that the luminosities inferred from the PC interpre- 
tation would be compatible with those predicted by equa- 
tion (3). In addition, as Arons (1981) emphasizes, equation 
(3) gives a minimum expected luminosity — although the 
luminosity of PSR J0437-4715 coincides with that mini- 
mum, PC luminosities of other pulsars may be greater. 
Thus, the soft X-ray radiation of at least some other mil- 
lisecond pidsars can be partly supplied by their PCs, al- 
though a nonthermal component may also be present and 
may even dominate. To resolve the dilemma, and to sepa- 
rate the two components, further observations arc needed. 
At present, we can only say that both possibilities should 
be explored, and complementary interpretations are still 
viable for some pulsars. 

6. Summary and conclusions 

We have shown that both the spectra and the light curves 
of the soft X-ray radiation of PSR J0437-4715 observed 
with the ROSAT PSPC and HRI and EUVE DSI can 
be interpreted as thermal radiation from two hydrogen- 



covered polar caps. The simplest, single-temperature PC 
model allows us to estimate typical PC radius, i?pc ^ 
1 km, and temperature, Tpe ~ 1 x 10^ K. The success- 
ful fitting of the data with the two-temperature model 
indicates that the PC temperature may be non-uniform, 
decreasing from Tcore ~ (LI — 1.8) x 10^ K at a PC core 
with a typical size i?eore ~ 150 — 400 m down to 
(3 — 5) X 10^ K at a much greater radius i?rim ^2 — 6 km, 
so that the heated area may comprise a considerable frac- 
tion of the NS surface. The bolometric luminosity of the 
two PCs is Lboi = (1.0 - 1.6) X lO^Orffgo erg s~\ that is 
~ (2 — 4) X 10"** of the total energy loss of the pulsar is 
absorbed and re-emitted by the PCs. 

We emphasize that the proposed interpretation is 
based essentially upon the neutron star atmosphere mod- 
els — the data cannot be interpreted as purely thermal 
PC radiation with the simplistic blackbody model. More- 
over, with the aid of the NS atmosphere models we show 
that the emitting layers should be depleted of heavy ele- 
ments, which can be naturally explained by the plausible 
assumption that the old NS has experienced accretion of 
the hydrogen-rich matter during its long life-time. 

The successful fits of the EUVE and ROSAT hght 
curves in different energy ranges with the PC models 
are possible only if the energy-dependent limb-darkening 
of the hydrogen atmospheres is taken into account. The 
shape of the light curves is also very sensitive to the bend- 
ing of photon trajectories in the strong gravitational field 
of the NS and to the orientations of the rotational and 
magnetic axes of the pulsar. We obtained the satisfactory 
fits assuming Mns = 1.4M0 and i?NS = 10 km, with the 
magnetic inclination angle, a = 35°, and viewing angle, 
^ = 40°, inferred from radio observations. We have demon- 
strated that the effect of the gravitational bending on the 
light curves can provide important information on the NS 
mass-to-radius ratio. 

The hydrogen column density towards the PSR J0437- 
4715, Uh ^ (1^3) X 10*^ cm~^, as well as strong ionization 
of hydrogen, ^ > 20%, estimated from the PC model fits 
are consistent with the ISM properties obtained from ob- 
servations of other objects in the vicinity of the pulsar, 
whereas the power-law fit yields greater uh and lower ^. 
This can be considered as one more argument in favor of 
thermal origin of the X-ray radiation of PSR J0437-4715. 

The inferred PC temperature, radius and, especially, 
luminosity are in excellent agreement with the predictions 
of the slot-gap pulsar model (Arons 1981). The upper 
limit on the PC luminosity provided by the outer-gap pul- 
sar model (Cheng et al. 1986) exceeds the observed value. 
This model becomes compatible with our results if a re- 
duced efficiency of the particle accelerator, consistent with 
the upper limit on the gamma-ray flux from the pulsar, is 
assumed. 

Radio pulsar models predict that heated PCs is a com- 
mon phenomenon inherent to all active pulsars. By virtue 
of this, the dependence L^; ^ O.OOIL; found by BT97 for 26 
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pulsars detected with the ROSAT determines a fraction 

of the rotational energy loss rc-cmittcd in X-rays by both 
nonthermal and thermal processes. Contributions of these 
two mechanisms to the total X-ray flux depend on vari- 
ous factors, such as pulsar period, radius, magnetic field, 
etc. In some pulsars (e. g., PSR B1821-24) the thermal- 
like soft X-ray radiation from PCs may happen to be less 
intensive than the nonthermal radiation generated in the 
pulsar magnetosphere. For instance, both the slot-gap and 
outer-gap models predict the luminosity of PCs of PSR 
B1821-24 much lower than the luminosity observed with 
the ROSAT and ASCA. 

Although all the available observations of PSR J0437- 
4715 arc consistent with the suggested PC interpretation, 
and this interpretation is supported by the theoretical pul- 
sar models and by indirect arguments, we still cannot 
completely exclude that the same data, collected in the 
relatively narrow energy range, could be interpreted with 
another model — for instance, as a combination of the 
PC radiation and magnetospheric radiation. Hence, the 
proposed interpretation can be considered as complemen- 
tary to the nonthermal interpretation discussed by BT97. 
A crucial test would be provided by observations of this 
object at energies > 1 — 2 keV, which would enable one to 
firmly discriminate between the power-law and thermal- 
like spectra, or to separate their contributions. Such obser- 
vations could be carried out by ASCA (with a sufficiently 
long exposure) and by the forthcoming AXAF, XMM 
and ASTRO-E missions. Important data could be also 
obtained from observations of the pulsar in the UV range 
(1200 < A < 3000 A) with the HST. In particular, such 
observations would enable one to measure the tempera- 
ture of the entire NS surface (cf. Pavlov, Stringfellow & 
Cordova 1996a) and to elucidate other possible heating 
mechanisms competing with the PC heating. 

For more accurate and reliable interpretation of the 
future observations, more theoretical work is highly desir- 
able. In particular, the temperature distribution over the 
NS surface within and around the PCs should be investi- 
gated (and further used to fit the data), and the nonther- 
mal X-ray spectrum and beam shape should be calculated. 
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